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Abstract 

A canonical approach to Cherenkov radiation in the presence of a magnetodielectric medium is 
presented in classical, nonrelativistic and relativistic quantum regimes. The equations of motion 
for the canonical variables are solved explicitly for both positive and negative times. Maxwell and 
related constitute equations are obtained. In the large-time limit, the vector potential operator 
is found and expressed in terms of the medium operators. The energy loss of a charged particle, 
emitted in the form of radiation, in finite temperature is calculated. A Dirac equation concerning 
the relativistic motion of the particle in presence of the magnetodielectric medium is derived and 
the relativistic Cherenkov radiation at zero and finite temperature is investigated. Finally, it is 
shown that the Cherenkov radiation in nonrelativistic and relativistic quantum regimes, unlike its 
classical counterpart, introduces automatically a cutoff for higher frequencies beyond which the 
power of radiation emission is zero. 
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I. INTRODUCTION 



Cherenkov radiation is the radiation with continuous spectrum and specific angular dis- 
tribution emitted by the medium due to the motion of charged particles moving in the 
medium with a velocity exceeding the phase velocity of light in the transparent medium. 
The Cherenkov radiation in transparent media was experimentally observed by Cherenkov 
1[. Theoretical explanation on this phenomenon was first developed by Frank and Tamm 
2). They showed that the particle should radiate when its velocity exceeds the velocity of 
light in the medium as the emitted rays make an angle 9 with the charge velocity given by 
cos# = c/vn where v is the speed of the particle, n is the index of refraction of the medium. 
Also they are polarized with the electric field in the plane of this angle. The radiation 
chock-front, called the Cherenkov cone, is analogous to the Mach cone formed as objects 
move with supersonic speeds through air. The cherenkov theory has drawn a great deal of 
attention all around the world. This theory has been widely used in high-energy particle 
physics, optics, cosmic-ray physics, h.gh-power Ration sources, aad so on gg Typ.ca. 
examples are the discoveries of the anti-proton [7] and the J-particle |8|. 

The classical theory of Cherenkov effect is sufficiently accurate in the optical part of spec- 
trum. For methodological reasons, it is equally important to consider the quantum theory 
of this effects. The phenomenological quantum theory of the Cherenkov radiation developed 
by Ginsburg but the dissipative character of the uredium was negated fl-Q. The source- 
theory explanation of this effect was given by Schwinger et al. [11] . Unfortunately, due 
to the method of combining the denominators of the propagators in parametric form, the 
resulting integrals are exceedingly complicated and approximations were necessarily made. 

The first formulation of the finite-temperature quantum field theory was presented by 
Dolan and Jackiw [l^ . Weinberg 13] and Bernard [ijj and the first application of it con- 
cerned the effective potentials in Higgs theories. Also, the inclusion of temperature has 
Deen carried out in QED of Cherenkov radiation only in a nondispersive dielectric medium 
I5I ]. The main purpose of the present work is to develop a canonical theory of the finite 
temperature Cherenkov radiation to evaluate the electromagnetic field arising from the uni- 
form motion of a charged particle in the presence of a magnetodielectric medium in different 
regimes i.e, classical, non-relativistic and relativistic quantum regimes. To achieve this goal 
we first generalize a Lagrangian introduced in [3] to include the external charges. This 
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prepares not only the grounds to extend the Ginzburg theory to include the dissipative 
character and the permeable character of a medium but extracts a Dirac equation for a rel- 
ativistic particle in the presence of the magnetodielectric medium to survey the relativistic 
effects of the Cherenkov radiation to finite temperature regime. 

The layout of the paper is as follows: In section 2, a Lagrangian for the total system is 
proposed and a classical treatment of the Cherenkov radiation to finite temperature regime is 
investigated. In sections 3, we use the Lagrangian introduced in the section 2 to canonically 
quantize the electromagnetic filed arises by the moving external charges embedded in the 
magnetodielectric medium. Maxwell and constitute equations are obtained. We find that 
for sufficiently large times the vector potential operator can be expressed in terms of the 
initial medium operators. The consistency of these solutions for the vector potential operator 
depend on the validity of certain velocity sum rules. We also show how to relate the results to 
the damping polarization model and phenomenological quantization theories. By considering 
finite temperature effects, energy loss of a charged particle emitted in the form of radiation is 
calculated. Subsequently, this formalism is generalized somehow to describe the relativistic 
moving particles. Finally, we discuss the main results and conclude in section 3. 



II. CLASSICAL THEORY 



Cherenkov radiation has the property that it occurs with uniform motion of a charged 
particle in a spatially homogeneous medium At first we attempt to treat the theory 
of Cherenkov radiation in the presence of a linear homogeneous magnetodielectric medium 
on basis of the classical theory. In the first part of this section we generalize the approach 
presented in |l6] to the case where there are some external charges in the medium and in the 
second part we examine the emission of electromagnetic wave, if it occurs, by the moving 
particle. In the third section we will concentrate our attention on the Cherenkov radiation 
in the finite temperature situation. 



A. Classical dynamics 

Classical and quantum electrodynamics in a linear magnetodielectric, can be accom- 
plished by modeling the medium with two independent reservoirs that interact with elec- 
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tromagnetic field. Each reservoir contains a continuum of three dimensional harmonic os- 



cillators describing the polarizability and magnetizability properties of the medium [18]. 
Therefore, in order to have a classical treatment of electrodynamics in a homogenous mag- 
netodielectric medium, we begin with the following classical Lagrangian for the total system 
(medium+electromagnetic field+external charges) 

L(t) = L res + L em + L q + Li n t. (1) 

The first term L res is the reservoir part 

Lres = J ^y"d 3 ri[X w (r,t)-X w (r,t)-^ 2 X w (r,t)-X w (r,t)] (2) 

+ J dcoj rf 3 r^[Y w (r,t)-Y w (r,t)-u; 2 Y w (r,t)-Y w (r,t)] 

where the dynamical variables X^(r, t) and Y^(r, t) correspond to the electric and magnetic 
characters of the medium, respectively. The second term L em is the electromagnetic field 

L em = /dM^oE 2 (r,t)-5!M]. (3) 
The third term L q is the Lagrangian of the external charges with mass m a and position r a 
L q = i^m a r 2 (t) + ^2(q a r a ■ A(r a ,t) - q a <f)(r a ,t)), (4) 

a a 

and finally L int is the interaction term, which includes the linear interaction between the 
medium and electromagnetic field trough coupling functions f{oo) and g{u>) and also the 
interaction between the external charges and electromagnetic field 

/OO f /*oo f 

dcu J d 3 r /(a;)X a) (r, t) ■ E(r, t) +J du J dh g(u)Y u (T, t) ■ B(r, t) (5) 

In the equations ([3]) and ([5]), E = — ^ — V0 and B = V x A are the total electric and 
magnetic fields and A and are the vector and the scaler potentials. For simplicity we work 
in the reciprocal space and write the fields in terms of their spatial Fourier transforms. The 
range of the variable k in the reciprocal space is restricted to the half space denoted by a 
prime over the integral i.e., J c? 3 k [19|], thus in the reciprocal half space the Lagrangian ([1]) 
can be written as 

L(t) = Ut) + L res (t) + L em (t) + L mt (t) (6) 
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L, = I E m ^ W + (Jp E / rf3k fc» (*« • ^( k ' *) - *)) elk rQ + c - c -] > (7) 



OO /■' /"OO 

3WIY 12 ,..2| Y |2\ , / J, , / j3WIv |2 , ,2| V |2N 



ciw / tfk(|Xj" -w 2 |Xj") + / I d 6 k(\Y 

ii y Jo 



ur Y, 



rf 3 k(e |A| 2 + e o |k0| 2 - l k X -I 2 ) + e / d 3 k(-?k • A0* + ex.; 



/'o 



(9) 



-it it 



doj I crk 



f(w) X£(k, t) • («Mk, t) + A(k, t)) + c.c. 

pOO n' 

+ J du J d 3 k [g(oj) 3C(k,t) • (zk x A(k,t)) + c.c] 



(10) 



where we have applied X*(k, t) = X(— k, i) and the similar relations for the other dynamical 



fields 



191 ] . we can obtain the classical equations of the motion simply from Euler-Lagrange 



equations. For the vector and scalar potentials A(k, t), 0(k, t) we have 



d 



5L 



5L 



i = 1,2,3 



dt \8(A*(k,t)) ) 5(A*(k,t)) 
=> yU eoA(k, t) + /i O e o «k0(k, i) - k x (k x A(k, £)) 

// P(k, t) + ^ k x M(k, t) + hq J(k, £) 



(11) 



and 



d 
eft 







k ,t))y W(kt)) 

-e «k ■ A(k, t) + e o k 2 0(k, i) = -zk • P(k, t) + p(k, *) 



where 



E(k, t) 



dojf(u)^(k,t), 



M(k,t) = / dug(u) Y^k,*), 
Jo 



J(M) 



p(k,t) 



-W2 E 9a *" Q 



(2tt) 3 / 2 



^3/2 E ?Qf 



-ikr Q 



(12) 

(13) 
(14) 

(15) 
(16) 



(2tt) 3 /2 

are respectively the Fourier transforms of the electric and magnetic polarization densities of 
the medium. The Fourier transforms of the external current and charge densities satisfy the 



charge conservation in the reciprocal space p = —tk ■ J. Similarly from the Euler-Lagrange 
equations for the fields r a and X^, we find 

d ( 8L \ 5L 



dt \S(r a>i (t))J 6(r a ,i(t)) ° 

m a T a (t) = q a E(r a ,t) + q a r a x B(r a ,t) 



i = 1,2,3 



(17) 



and 



5L 



5L 



dt W(xl(k,t)W <*Q&(M)) 



i = 1,2,3 



X^(k, t) + w 2 Xjk, t) = -f(u) (A(k, t) + ik<p(k, t)) 



(18) 



5L 



5L 



dt UC&fct))! SCCiiKt)) 



% = 1,2,3 



X,(k, t) + w 2 Yjk, t) = g{u) ik x (A(k, t). 



The formal solution of the field equation (|T8|) is 

/, x ,v /, ^sincot „, , f ,.smcj(t — t') 

20k,i) = X (k,0) + X (k,0 coso;t + / w / di' * ^ 

— u - Jo 



Efk,n 



(19) 



(20) 



where the first term is the inhomogeneous solution of the equation ({TBI and the second term 
is the homogeneous one. We will show that after quantization the homogeneous solution 
becomes a noise operator. However, since we are only interested in the induced polariza- 
tion, we keep the inhomogeneous solution and using the equation (fT3"|) . we find the electric 
polarization density of the medium in reciprocal space 



/"OO 

P(k,t) = e / d*Xe(*-OE(k,0 
Jo 



(21) 



where Xe is the electric causal susceptibility of the medium and in terms of the coupling 
function / can be written as 



J, r^ E^ fy 

eo JO u <> \ / 



t > if 



Xe(t -t')= I 



(22) 



t < if 
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which is the origin of the significant Kramers- Kronig relations [17|]. In a similar fashion the 
magnetic polarization density of the medium in reciprocal space can be obtained straight- 
forwardly using the formal solution of the equation (1T9"1) as 

1 f°° 

M(k, t) = — dt Xm (t - OB(k, (23) 
A^o Jo 

where \m is the magnetic causal susceptibility of the medium which in terms of the coupling 
function g can be written as 

\ Q ^dw^^lg\u) t>t> 

(24) 

t < t' 

The electric permittivity and the inverse magnetic permeability of the magnetodielectric 
medium are defined in terms of Xe and Xm as 



Xm(t-t') 



(w) = l + Xe(w) (25) 



and 



where 



k(u) = 1 - XmH (26) 



DO 



JO 

By using Eqs. (|22|) and ff24l) . we can obtain the following important relations in the frequency 
domain 

Xe(co) = - f°°dw' a /2( f } (27) 

/"°° o 2 fw / ) 

Xm{co) = /io / ^ , 2 V n+ - ^ 28 ) 
Jo u/ 2 — w 2 — z0+ 

These are complex functions of frequency which satisfy Kramers-Kronig relations and have 
the properties of response functions i.e, e(— w*) = e*(co) and k(— w*) = n*(cu) and Ime{u) > 
0, Imfi(u) > provided that f 2 (—u*) = f 2 {oo) and g 2 (— u*) = g 2 (u). It can be shown that 
these functions have no poles in the upper half plane and tend to zero as u — > oo. If we are 
given definite electric permittivity and inverse magnetic permeability of the medium then 
we can inverse the relations ( |22l) and (1241) and find the corresponding coupling functions 
f{uj) and g{uj) as 



fU)= \/^Ime(u) (29) 



g(cu) = \ IrriK(uj) (30) 

V 

where the minus sign in the second expression is used since for a magnetodielectric 
medium Imfi(u) > 0, therefore Iitlk{uj) < 0. In order to illustrate the relations between 
the coupling functions and the electric permittivity and the inverse magnetic permeability 
of a magnetodielectric medium, let us restrict our attention to a single resonance electric 
permittivity which can be obtained from the Lorentz oscillator model 

= 1 + t— ^ (31) 



and a single resonance magnetic permeability [20] 

M") = i + - — If — (32) 

where u pe and u pm are the coupling strengths, uo m , uo m are the transverse resonance fre- 
quencies, and 7 e , 7 m are the absorption parameters. Now by using Eqs. (|29|) and (130]) the 
coupling functions f(u) and g(u) can be obtained as 

fit \ 2 le e u 2 u 2 /n 

(^0e - W 2 ) 2 + 7 2 ^ 2 

2/ \ 2 7m ^ m q; 2 /7r/io 

# M = 7~~2 i 2 2 N 2 , 2 2 - ( 34 ) 

Following the standard approach to classical electrodynamics, we choose the Coulomb 
gauge k ■ A = 0, so that the vector potential A is a purely transverse field. By useing the 
Euler-Lagrange equation for if*, we eliminate ip from equation fill I) and substituting Eqs. 
fl2"T|) and fl2"3"|) into Eq. (fTTj) . obtain the following inhomogeneous wave equation 

f* d f* 

^ eoA(k, t) + /c 2 A(k, t) - m>k 2 / dt' Xm (t - t') A(k, t') + fi Q — / dt^t - A(k, t') 

= /i J ± (k,t). (35) 

where the transverse current is defined as J -1 = Y^x=i^-' e A(k) with unit polarization vectors 
e^(k),A = 1,2, which are orthogonal to e^(k) = jr = k and to one another. This equation 
can be solved in terms of initial conditions using the Laplace transforms. For any time- 
dependent operator Q(t) the forward Laplace transform is defined as 

POO 

tt f (s) = / dte- st tt(t), (36) 
Jo 
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obviously, the Laplace transform contains all information about the time evolution of Q for 
positive t. In the following we wish to determine the time evolution of the relevant operators 
of our model for any time, either positive or negative. Hence, we also introduce the backward 
Laplace transform 



oo 



n\s) = / dte- st Q(-t). (37) 
Jo 

Let e(s) and k(s) be the Laplace transformations of e(t) and K{t), respectively. Then A (k, s) 
and A (k, s) i.e., the forward and backward Laplace transformation of A(k,t), can be ob- 
tained as follows 

~ (k ' S) = s 2 ~e(s) + k 2 c 2 k(s) MK 0) ± A-( s ) + W(i) A(k ' 0) 
1 v (J / - b (k, g )-e A (k))e A (k) 

e V s 2 e(s) + k 2 c 2 k(s) 1 J 

The time-dependent vector potential is obtained from a contour integration over the 
Bromwich contour by an inverse Laplace transformation. From Eq. (138 j) we find the vector 
potential for t > 0, [2lJ 

A(k,t) = e(*U(k,0) + C(t)i(k,0) 

2 — E r^ a f (^;»H>w) ei(k) , (39) 

27re ^ J-oo w^e(w) - k 2 c 2 K(u) 



where 



and 



1 p-\-tOO ~i \ J 

= -L / dse st S£ i s = y ^(e-rfij*^) (40 ) 



C(*) = — / ds^-r— ——, — — = — Ylmie-^— |— ) (41) 
sw 2mJ_ loo s 2 e{s) + k 2 c 2 k(s) kc^ v cfv(^) ; v ; 

In these expressions, we changed the integration variable from s to — iuj + 77, with a small 
but positive 77. Therefore, we introduce the electric permittivity and the inverse magnetic 
permeability of the magnetodielectric medium in the frequency domain as e(u) = e(— tu + 0) 
and k(u) = k(—iu + 0) for real u, and correspondingly, e(tu + 0) and k(iu + 0) as their 
complex conjugations e*(u) and k*(u), respectively. We define for each allowed frequency 
Qj, the group velocity v J g = || and the phase velocity v J p = ^ where the frequencies %(k) 
and fi?(k) are the complex-frequency solutions of the dispersion relation w 2 e(u;) — k 2 c 2 /t(u;) 
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which has no zeros in the upper half plane. It is worth emphasizing that for a lossy medium, 
we lose the usual dispersion relation in which a limited number of discrete frequencies u are 
associated with each wave vector k. Thus k and u must be considered as independent real 



variables 



22] 



Now the vector potential for t < is obtained from the inverse Laplace transform of Eq. 
fl38]) as 



A(k,t) = £(*)A(k,0)-C(t)A(k,0) 

zTieo ~^ J-oo u 2 e*{u) — k z c z K*(u) 

A 

The coefficients rj(t) and in Eqs. and ( I42p damp out exponentially in time since 
all Oj in the exponentials have negative imaginary parts. Therefore, For large times, the 
medium and electromagnetic field tend to an equilibrium which is determined by the charac- 
teristic damping time Tj = 1/ImQj. After a few times the maximum characteristic damping 
time, only the third term survives in these equations since they have poles on the imaginary 
axis in the complex s plane. Also, these terms in Eqs. (139]) and (142]) are zero for negative 



and positive t, respectively 



231 ] . Thus, we may combine the two expressions into a single 



one and use the Fourier transform to obtain the vector potential in real space for all t 

A(r,t) = y du J d 3 ke ik - r -'^A + (k,w) + c.c., (43) 
with the positive-frequency Fourier component 

A + (k.^) = [ - ( ,Yr + °V f\ ] + fe/TS \ e t\ wv (44) 

(27r) t, / z e ' w z e(u;) — k z c z n(uj) u) z €*(uj) — k z c z K*{u) 

A 

It is not difficult to show that the transverse electric field can be written as 

/+oo /• 
duu J d 3 k{e lk - r - lult A + {k,uj) - c.c.) (45) 

From the Eq. (139|) we can see that the coefficient £ (t) in Eq. ([40]) takes the value 1 at time 
t — and therefore the coefficient takes the value 0. These constraints are satisfied, if 
certain velocity sum rules are adopted. In this way, we find the following modified velocity 
sum rules for all wave vectors k 

£ife(4) = l (46) 
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and 



0. 



which are resembling the quantum relations obtained in 
magnetodielectric medium. 



(47) 

2l| and 24j, now generalized to a 



B. Classical theory of Cherenkov radiation (T=0) 

Theoretically, when considering the Cherenkov radiation, one usually treats the charge 
motion with a constant velocity which corresponds to the so-called Tamm Frank problem [2j . 
Consider a point charge e uniformly moving in a magnetodielectric medium with a velocity 
v. Therefore, according to Eq. (jI5T) 

P V 

I(M)= (2^p e ' (48) 



then 



j> (k , -«„ + o) = * ( w - k • v ) • ( 49 ) 

Now substituting Eq. (|49|) into the equation (j4"5j) . we obtain 

E ± ( r , t) = V ! +0 ° douu [ d 3 k . e^'-e ^v-e A (k) _ _ 

v ' ; 87T 3 e V 7 [ w 2 e(w)-Fc 2 K(u) J v ; v ; v ; 

A 

In this case the energy loss of a point charged particle per unit length emitted in the form 
of radiation, is defined by the braking force acting on the charge at its location {25 ] 

dW 



ev ■ E 

le 



r=vt 



dt 

r+00 , .(k-v-w)*^ . eA ( k ))2 



(51) 



87r 3 e J J k 2 c 2 K,(u) — u 2 e(cu) 

A 



Letting 6 be the angle between v and k, then X1a( v ' e A(k)) 2 = v 2 (l — cos 2 6). Therefore we 
find 

dW ie 2 v f + °° , f + °° „, f +1 u XkvcosO-^tn _ cos 2m w 

^ = ^%X ^y dfcfc I, rf(cosg)[ - ~ c - c w°° ag ~ ^ 

e 2 ^ f + °° , /"+ 00 „ , Im(k 2 c 2 K(uj) - w 2 e(w)) . w 2 . 

' duu, dkk A > ) >\ i- ) (52) 



27r 2 e Jo yo |/c 2 c 2 k(o;) — u; 2 e(ci;)| 2 /c 2 v 2- 
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that the electromagnetic waves are emitted at an angle to the path of the particle determined 
by 

cos6 = (53) 
kv 

The transparent magnetodielectric medium can be considered in principle as a limiting case 
of the lossy dispersive medium. Actually, there are some ranges of frequencies over which 
the imaginary parts of the permittivity and permeability of the magnetodielectric medium 
can be ignored. In a transparent magnetodielectric medium the energy loss of a charge a 
only a result of radiation. In this case the imaginary parts of the electric permittivity and 
the inverse magnetic permeability must tend to zero, thus 

Ime{u), Jm/x(w)-»-0 \k 2 C 2 K{UJ) - L0 2 €{UJ) \ 2 



(54) 



where n 2 (u) = e(u)fi(u) and the frequencies fij(k) are the complex-frequency solutions of 
the dispersion relation u 2 e(u) — k 2 c 2 K(uj). Now Substituting Eq. ([51]) into Eq. (|52|) and 
performing the integration over u and converting the integration over k to an integration 



over Qj, we obtain 



r\W p 2 d -, f + °° r 2 

E / ^'M^)(1--™t) (55) 



It is seen from the Dirac 5-function in Eqs. (154]) and (]52]) that radiation is possible only 
if the inequality v > c/n(Qj) is satisfied. We define the Cherenkov cone cosO = c/vn(Qj) 
corresponding to any frequency Qj for which v > c/n(Qj), where 6 is the angle between the 
wave vector k of the radiated electromagnetic wave and the velocity of the particle v. It is 
easy to show that when there is only one frequency for each k, then the equation ( 155]) tends 



to the result of 



9|, (l7| and |26|. 



C. Finite temperature Cherenkov radiation in classical regime 

Our considerations so far have been applied to zero temperature. The generalization of the 
formalism for this case is straightforward. It is known that the medium and electromagnetic 
field are in the thermal equilibrium in this regime. The inclusion of temperature may be 
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done in the usual manner [27]- [29]. The finite temperature expression, as is well known, is 
found by replacing the frequency integral by a sum over Matsubara frequencies according 
to the transition 

^-Mi) k B Tj2f«i), £i = 2nk B Tl/h (56) 

n 1=0 

where T and k B are the temperature and Boltzmann constant and the prime on the sum- 
mation mark denotes that the zeroth term is given half weight as is conventional. The effect 
of finite temperature on the energy loss in the form of Cherenkov radiation can be easily 
taken into account by using Eqs. (156]) and fl52|) 

dW ie 2 vk B T ^ f r^^ ImiWKte) J 2 
dt ne h 2^*J \k*<*K(&) + $e(tb)\* 1 + kW h { *' } 

We know that the function coth(hu /2k B T) has an infinite number of poles at uji = i& and 
elsewhere is analytic and bounded. This enables us to write Eq. (1571) as 

dW e 2 v f + °° , w Huj . f +QO „ Tm{k 2 c 2 K(uj) - u 2 e(u)) „ w 2 , ro , 

Then with suitable rearrangement of the exponentials in the hyperbolic cotangent we obtain 

dW f dW s ,dW. 

)t=o + (-jtWo ( 59 ) 



where 



(it rft (it 



e 2 w /-+ 00 /• +0 ° Im(k 2 c 2 K(u) - u 2 e(u)) „ w 2 



dt T 27r 2 eo Jo */o |/c 2 c 2 k(u;) — w 2 e(u;)| 2 v A; 2 f 2 



dkk ^T^Z77A ..»./..m, C 1 -?^) ( 6 °) 



and 



(W e 2 ?; /- +00 2a; / + °° wt-t. Im (k 2 c 2 n(u) - m 2 e(tu)) u 2 

{ ~aT }T ^~2~^ 2 ToJo e^l 2k BT)-i} \k 2 c 2 K(u)-u 2 e(u)\ 2 ( 

(61) 

The last formula differs from the zero-temperature formula only by the multiplicative factor 
2[e <yhw l 2kBTS) — 1] _1 which has the asymptotic behavior and Ak B T/t\oj at low temperatures 
k B T « fajj and at high temperatures k B T » hu, respectively. In fact, the Matsubara 
frequency sum naturally separates into a term which is temperature independent and a term 
containing the Bose-Einstein distribution. In some sense the replacement of Matsubara 
frequency sum by an integral, as in (I58|) . is equivalent to switching from imaginary time 
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(discrete frequencies in Euclidean space) to real time (continuous energies in Minkowski 
space) . 

A transparent magnetodielectric medium can be considered in principle as a limiting case 
of a lossy dispersive medium. In this case using Eq. (|54p . we obtain 

dW e 2 v f + °° „ , , ML c 2 , 

which is the finite temperature generalization of Eq. (j55l). It is easy to show that for a 
nondispersive medium, Eq. (J55l) tends to the result of |3Ql ] - 

III. QUANTUM THEORY 

The classical theory of cherenkov radiation effects is sufficiently accurate in the optical 



part of the spectrum lOj. For methodological and physical reasons, it is equally important 
to consider the quantum theory of this effects. Quantum theory enables us to derive the 
classical equation with the appropriate corrections. We extract the Maxwell equations and 
constitute relations and the vector potential field operator in the first part of this section, and 
in the following we calculate the radiation intensity within a nonrelativistic and relativistic 
theory. 

A. Canonical quantization 

In the description of the canonical quantization of the electromagnetic field, we choose the 
Coulomb gauge k ■ A(k, t) = 0. In this gauge the vector potential A is a purely transverse 
field and can be decomposed along the unit polarization vectors e,\(k) A = 1, 2 

2 

A(k,t) = J]A A (k,t)e A (k). (63) 

A=l 

The dynamical fields an d has both transverse and longitudinal parts and can be 
expanded as 

3 

2UM) = ^2LA(M)e A (k). (64) 

A=l 
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Furthermore, by using the Lagrange's equation for the scalar potential 0, we find a La- 
grangian depending on a reduced number of dynamical variables in the reciprocal space 



a 

+ Hj duj J d3k (\X,x\ 2 - ^ 2 \2Lx\ 2 + \t,x\ 2 - ^ 2 \y, 



\2\ 



Al 



A=l 
2 



/\kA r f 
rf 3 k(e |i A | 2 - ±=*-) + / rf 3 k(A A Jf (k,t) + /i.e.) 
/io 7 

+ ^ / C / 3 k[(-i A P A ± (k,t) + (zkx ( i 4 A e A (k))-M*(k,t) + / i .c] 

A=l 





pi 


2 




|A 
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+ / d 3 k(- + h.c.) - / d"k^=— - / rf 3 k ' .-' . (65) 

The Lagrangian fl66|) can now be used to obtain the corresponding canonical conjugate 
variables of the fields within the half k-space as 

-RUKt) = = e A x (k,t) - Pt(k,t) (66) 

R^ k ' *) = TT^T = ^a, ILA(k, t) = = t- a (67) 

Pa(*) = tt^t = "\*r Q + g Q A(r a , £). (68) 

(?a ) 

Now the fields can be quantized canonically in a standard fashion by demanding equal-time 
commutation relations among the variables and their conjugates. For electromagnetic field 
components, and the dynamical variables of the external charges, we find respectively 

[Al{k,t), -D_i,(k',t)] = ihS xy 5{k - k') (69) 

[r a i(t),p aj (t)} = ihSij (70) 

and for the reservoir fields 

x: A (k,t),Q^, A ,(k',t)] =ih6 xy 6(uj-u')6(k-k') (71) 

[F: A (k,t),rL, v (k',t)] =ih5 xx ,5(cu-u')5(k-k') (72) 

with all other equal-time commutators being zero. Using the Lagrangian (|66|) and the 
expression for the canonical conjugate variables in fl69|) . we obtain the Hamiltonian of the 
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total system 



A=l 

3 coo 



2m a 



+ £ fW rf3k (il.Ai 2 +^i^Ai 2 + ^Ai 2 + ^ir.Ai 2 ) 

A=l ^ J 

- X) / rf3k t( zk x UA e A(k)) -M*(M) + ^.c] + / rf 3 k 

A=l ^ ^ 



zk ■ PI 2 



eol^l 2 

. , . -ik-P)p* f _ |p| 2 , 



e \k\ 

If we apply Heisenberg equation to the operators D_ x and A x , and use the commutation 
relation (EH]), Maxwell equations in the reciprocal space can be obtained as 

= i[fr,4x(M)] = - - f(k,t) " £f(k,t) (74) 

Ai e 

Z^(k,t) = i[if,Dj!-(k,*)] = — 4(k,f)-e A (k).( ? kxM(k,i))- J^(M) (75) 
ft po 

Multiplying both sides of these equations by the polarization unit vectors and summing over 
the polarization indices we find 

B X (k, t) = eoE^ (k, t) + P x (k, t) (76) 



D (k, t) = ik x H(k, t) - J 1 (k, t) (77) 

where D -1 is the transverse displacement field, E -1 = —A is the transverse electric field and 
/xoH(k, t) = ?k x A(k, i) is the magnetic induction field and 

J X (k, t) = E + e~^r a ). (78) 

In the presence of external charges, the longitudinal components of the electric and displace- 
ment fields can be written respectively as 

D" (k, t) = e E" (k, t) + P" (k, t) = ^fi . (80) 
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If we differentiate the equation (1761) with respect to the time variable t and use the equation 
( 177j) . we find the quantum counterpart of the equation f fTTl) for the vector potential in the 
reciprocal space as 

/i e A(k, t) + |fc| 2 A(k, t) - At P ± (k, t) - n ik x M(k, t) = ^ Q J J -(k, t) (81) 

Equation flHTj) is the Langevin equation for the vector potential A(k, t), wherein, the explicit 
form of the electric and magnetic polarization densities of the medium is known. The 
quantum Langevin equation can be considered as the basis of the macroscopic description 



of a quantum particle coupled to an environment or a heat bath 3JJ. Similarly, it is easy to 



show that the Heisenberg equation of motion for the external charged particles is 

% 1 
m a r a = T[ H ,Pa -q a A(T a ,t)] = q a E(r a ,t) + -q a (r a x B(r a ,t) -B(r«,t) x r a ). (82) 



Using the commutation relations ([71]) . ([72]) and applying the total Hamiltonian (]73|) . it 
can be shown that the combination of the Heisenberg equations of the canonical variables 
X(u;,i) and Y(u,t) lead to the same equations (1181) and ( fT9l) with the solutions (1201) and 
( |23l) . respectively Now by Substituting the equations ( J20|) and ( 1231) in the integrands of the 
equations ( JT3|) . we find the the electric and magnetic polarization densities of the medium 
in reciprocal space 

POO 

P(k, t) = eo dt X e{t - t')E(k, + P iV (k, t) (83) 
1 F°° 

M(k, t) = — dt Xm {t - t')B_(k, f ) + M N (k, t) (84) 
Jo 

where Xe and x™ are the same electric and magnetic susceptibilities of the medium defined 
in the equations (T22l) and (T24l) . the noises 

roc • j. 

P iV (k,t)= / dw/(w)(2L(k,0) +2L( k >0)cosu;t), (85) 

M Af (k,t)= r dug (w)(Y w (k,0)^^ + Y w (k,0) cos ut) (86) 
Jo " w 

are the electric and magnetic polarization noise densities associated with absorption, with the 
causal behavior of the medium, respectively. To facilitate the calculations, let us introduce 
the following annihilation operators 



aA(M) = ^2^oWl ( e ° C ^l4A(k,t) -tDifrt)) , (87) 

17 



d x (k, oj,t) = J^- (ojX^ x (k, t) + *g wA (k, t) 



(88) 



b x (k,u,t) 



2t)w 



(u^ A (k,t)+4L A (k,t)). 



(89) 



From equal-time commutation relations for the fields (!69|) - (!72|) . we obtain the following 
equal-time commutation relations for the creation and annihilation operators 



a x {k,t),a[,{k',t) =5 xy 5{k-k' 



(90) 



d x (k, to, t), d x ,(k' , to', t) = 5 XX /5(co — u')5(k — k') 



(91) 



6 A (k, oj, t), 6 A ,(k', u', t)] = S xx ,6(u - u f )5(k - k') 



(92) 



The commutation relations fl90|) - fl92|) in contrast to the previous relations f l69|) -f J72|) . which 
were correct only in the half k-space, are now valid in the whole reciprocal space. Inverting 
the equations flH71) - fl88l) . we can write the canonical variables A, and in terms of the 
creation and annihilation operators as 

— 2 



A(k,t) 



h 



2e c\k\ ^ 



a A (k,t) + a A (-k,i))e A (k) 



X 



k, t) = J— ( rf A(k, u>, t) + 4(-k, u, t)) e A (k), 

U A=l 



YJM) = J]-Y1 (h(k,u,t) + b{(-k,uj,t^ e A (k) 

U A=l 



(93) 



(94) 



(95) 



Now by employing the Fourier transforms of these recent relations, the Hamiltonian of the 
total system (1731) . in the real space, can be recast into the final form 



H 



d s r\ 



D ± (r, t) ■ P(r, t) P 2 (r,t) 



fo 



+ 



2er 



2m n 



V x A(r,t)-M(r,t)] 

,(V-P(r,t))(V'-P(r',t)) 



r — r 



i3r (V-P(r,t)) | 1 ^ q a q p 



+ Hp + H e + H v 



r — r r 



+ P 



(96) 
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where 



P(r,t) = Y^J°° duJ j d3k \[^ f ^ K(k,c^)e lk - r + /ix.)e A (k), (97) 

A— 1 ^ 

M(r, t) = J2l°° duj J d " k (b x (k, w, i)c^ r + h.c.) e A (k) (98) 



and 



H F = J2[ d 3 k hc\k\ 4(k, £)a A (k, t) (99) 

A=l ^ 

H e = Y^ [ duj [ d 3 khw d{(k,uj,t)d x (k,u,t) (100) 

A=l ^ ^ 

# m = ^ /da; f d 3 khw 6 A (k, w, t)6 A (k, w, t) (101) 

A=l ^ ^ 

are the Hamiltonian of the electromagnetic field and the medium in the normal ordering 
form. 

We now proceed to solve the wave equation flHTl) along the lines of the classical wave 
equation (1551) using the Laplace transform technique. After some lengthy and elaborated cal- 
culations (see appendix), the vector potential in the large-time limit, i.e. when the medium 
and electromagnetic field tend to an equilibrium state, can be obtained as 



AM) = Z! y [ rf 3 k f ^ *_ m ( d ff[° )e ^'\ - h.c) e A (k) 



A=l 
2 



1 ^ f J3i fj / h k 2 / \ A & A (k, W, 0)c 6 , \ 

H > / cf k / duj , ' q(cu) \ ' oio , . - h.c. s A (k) 

eo^i 7 y2(27r) 3 u; yi ; \-co 2 e(tu) + c 2 k 2 n(oj) J XK) 

(102) 



where s A (k) = k x e A (k). In the large-time limit, the vector potential operator will be a 
function of the medium operators only, i.e, the radiation is due to the medium, which still 
satisfies Maxwells equations, as expected. Also the canonical commutation relation (I69p is 
preserved in this limit if in addition to the velocity sum rules (H6|) and ( H7|) . which still 
legitimate in the quantum domain, the following velocity sum rule for a magnetodielectric 
medium is also satisfied (see appendix) 

" +, ' x ui 3 Ime(ui) — k 2 c 2 ujImK{uj) n 



r+oo 

[A*(k,t),-D^(k',t)] = / du- 
L J J \— UJ 2 e(uj) + c 2 k 2 K(u)t\ 2 2 

^E^[^] = 1 - ( 103 ) 
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The form of vector potential operator given in fll02j) agrees with previous works, if one 
replaces the medium annihilation operators d\(k, u, 0) and b\(k, u,0) in the large-time fl96|) 
with the diagonalizing annihilation operators K ej \(k, u) and K mj \(k, co), derived by the 



damped polarization formalism 



24|, 32j and also if one makes similar replacements for the 



creation operators /|(k, u>) and /™(k, u) derived by the phenomenological formalism 33], 



34j . where again the same expressions for the field and medium operators are recovered. 



B. Nonrelativistic quantum theory of Cherenkov radiation 

We consider a charge particle with mass m and electric charge e uniformly moving in 
a linear homogeneous magnetodielectric medium described by the Hamiltonian (}96l) . In 
fact, we consider a total system of two noninteracting parts, that is, the free electron and 
a system, which consist of the electromagnetic field and the magnetodielectric medium in 
interaction. Therefore the Hamiltonian operator of the total system (196]) i.e electromagnetic 
field, the medium and the particle, in the large-time limit, can be rewritten as 



H = H + 



int 



Ho = H e i e + H F 

ele 2m 
3 



Hp-.y^l dco d 3 k frbj {d{(k,cu,t)dx(k,u,t) + b{(k,u,t)b\(k,u,t)} : 

A=l J J 

Hmt = _£>L^i) (104 ) 

m 

where x is the position operator of the particle and we have have ignored the term A 2 since 
the Cherenkov radiation can be considered as a first order process in which the number 
of photons changes by ±1. Also the direct Coulomb interaction between the electron and 
the medium has been omitted from the Hamiltonian (|96l) since it can give rise to radiative 
transitions only in third or higher orders. 

The moving charged particle with the momentum ^.q has the quantum state |q) 

p|q> = fiq|q>, (105) 

where |q) is the momentum eigenvector of the particle which in a coordinate representation 
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can be written as 

< x lq> = M*) = j^r^ ( 106 ) 

therefore the Hamiltonian H e i e = p 2 /2m for a free particle has the eigenvector |q) with the 
energy eigenvalues E q = K 2 ^ 2 /2m. The unperturbed Hamiltonian Hq = H e i e + Hp has the 
eigenstate 

| ele + rad) =\ ele) \ rad) 

which are the direct product of the eigenstates of H e i e and Hp. In order to separate out the 
emission of Cherenkov radiation from various other processes which might occur, such as, 
for example, ionization, emission of bremsstrahlung, etc., we restrict our attention to first 
order transitions. Therefore We apply quantum mechanical perturbation theory up to the 
first order approximation to treat the transition probability per unit time for a free particle 
of momentum ftq to emit a photon of momentum hk and energy huj thereby changing its 
momentum to /i(q — k) as following 

IW k = y |(l k I (q - k | H int | q) | 0) | 2 6(^£ -^\q-k\ 2 -hu) (107) 

where the states | 0) and | lk) present the vacuum state of the electromagnetic field and 
the excited state of the electromagnetic field with a single photon with wave vector k and 
frequency u, respectively. The argument of the Dirac 5 function displays the conservation 
of energy and the square of the matrix element of the equation (110 7ft is obtained by using 
Eqs. (PIED and ffTMj) as 

|(l k |(q-k|iJ m ,|q)|0)| 2 = - {|/3(a;,|fc|)| 2 |(q-k| e - k - x p-e A (k)| q )| 2 

+ fc 2 | 7 ( W ,|A;|)| 2 |(q-k| e - lk - x p-s A (k)|q)| 2 } (108) 

by using Eqs. (11061) and f ll(J5j) the matrix elements in above equation are just Kk ■ e A (k) 
and Kk ■ SA(k), respectively. Let 9 be the angle between q and k and let v = hq/m be the 
particle velocity, we find 



e 2 v(l — cos 2 #) { uj 2 Ime(u) — k 2 c 2 ImK{uj)\ „ uj . Hk 2 N 

\ 5 ^, — " — V 5 (cos 6 - — 1 + 

ATt 3 e hk V \-u 2 e(u) + k 2 c 2 K(uj)\ 2 J kv K 2mu' 



therefore, photon is emitted at an angle to the path of the particle given by 

uj , hk 2 , 

COS = — (1 + . 110 

kv K 2mu J K J 
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If the energy of the photon hw is much less than the rest mass of the particle mc 2 then 
this is approximately Eq. (l53j) which gives the classical Cherenkov angle. The total energy 
radiated per unit time is found to be 

dW e 2 v A 



_2„, 2 r r+°° 



dt ' 27r 2 e f- 



f i f 

/ d k / <ia;/ia;r q ^ q _ k 



A=l 

/ kdk / udu\l - —z — t (1 H rUm ^— = = ; : 



2n 2 e 



111) 



We note that the integration on the azimuthal angle is trivial. The integration on polar 
angle is done with the help of the Dirac S function in Eq. fillip and 

VM,)-^M,) \ = Im / i ") . (112) 

|-a; 2 e(o;) + A; 2 c 2 K(a;)| 2 / V-^M + k 2 c 2 K(u)J 



It is easily shown that Eq. fillip is consistent with the result of 
limit as H — > 0, reduces to the classical results ( !52l) . 



3511 and in the classical 



C. Relativistic quantum theory of Cherenkov radiation 



The description of particles used in the preceding section is valid only when the particles 
are moving at velocities small compared to the velocity of light. The preceding formalism 
must be generalized somehow to describe the relativistic moving particles, that is, we should 
drive the Dirac equation for external particles embedded in the magnetodielectric medium. 
For this purpose, we substitute the following Lagrangian for external particle instead of the 
Lagrangian (J3J) jsfj, {37] 



ihc 



W^(x,t) 7 A 

,u=0 ^ 



dx^ 



dx^ 



7^(x,t) 



— rac 2 ^)^) 



+ e / d 3 x 



t)7^(x, t) A J '(x, t)) - ^(x, t) 7 V(x, t)<p{x, t) 

.3=1 



(113) 



where 7^,/i = 0, ...,3 are the Dirac matrices with 7 = (3, 7 J = (3atj and ip = ip^(3. In a 
standard representation we have 



I > 
















[0-11 






/ 



(114) 
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where aj,j = 1, 2, 3 are Pauli spin matrices and I is the unit matrix. It is worth mentioning 
that the Lagrangian of the magnetodielectric medium fl8]) need not to be written in a covari- 
ant form since the medium is at rest and a non relativistic description is enough although 
it can be written in a covariant form straightforwardly. We proceed along the lines of the 
preceding section and instead of the canonical momentum of the particle p a we define the 
canonical conjugate variable of the Dirac particle ih^y as 

^ = (us) 

2 K ' 

The quantization procedure for the Dirac field can be achieved by imposing equal-time 
anticommutation relation among the field components 

{Va(x,t),^(x',t)} = 5 Q/3 5(x-x') (116) 

together with {^(x, t), ipp{x!, t)} = 0. Here we have chosen the anticommutation relation, 
since we are developing a theory of particles that obey Fermi-Dirac statistics. Using the 
Lagrangian ( 1651) and (11 13ft and the expression for the canonical conjugate variables in (1691) - 
(172"]) and (11161) . we obtain the Hamiltonian of the total system 

ff ^/A(-*,.(x, i )o.V.Wx, i ) +m c V (x, i )Wx, t) ) 
+ V / duo d 3 k( 

A=l J ° J 



—uj\ 



2 + u 2 |X lX | 2 + 



HwA 



2 



j d 3 k[{ik x Q4 A e A (k)) ■ M*(k,t) + h.c] + I d 3 k 

A=l •* J 



tk-P\ 2 



eo\k\ 



2 



, _ (-tk-P)p* f „ \p 

+ / d 3 k =^ + /i. c .) + / d 3 k - 



e |/c| 2 7 e |/c| 2 



- / d 3 k(J*(k,t)-A(k,t) + /i.c.) + ^ / ^( 9 ■ (H7) 

^ A=l ^ 2 

It is easily shown that the Heisenberg equation for the dynamic variable of the electromag- 
netic and medium fields lead to the same constitute equations (I83p - (|M|) and also Maxwell 
equations fl76l) - f!79l) in the reciprocal space. The Fourier transforms of the external current 
and charge densities are defined by 

J(k,t) = J^J~ 2 f rf 3 X^( X)t)aV ;( X)t ) e -^x (ng) 

p(M) = / rf 3 x^(x,t)^(x,t)e- 4k - x (119) 
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If we apply Heisenberg equation to the Dirac fields ^(x, t) and make use of the anticom- 
mutation relation (j!16p . Dirac equation in the presence of the electromagnetic field can be 
obtained as 



ihip(-x, t) = -cat ■ (iW + eA(x, t))^(x, t) + ey2(x, t) + mc 2 fi^{y., t), (120) 

where 



<P<* = 7^ ^ ( 7 1 ' + )e lk x , (121) 



-ik-P(k,i) p(M). 
(2tt) 3 / 2 J + €o k 2 

is the scalar potential of the electromagnetic field defined in Eq. ( fl2l) . We now expand the 
Dirac field ip(x.,t) in eigenfunctions of the Dirac equation in the absence of the electromag- 
netic field 

i 4 r 

^ (x,t) = j^fr^J rf VM(q^(q) (122) 

where ^(q) = tt M (q)e iq " x and w M (q) are four-component spinors of the Dirac equation with 
eigenvalues _E q = ±A/^ 2 c 2 q 2 + m 2 c 4 with the normalization condition (q)w„(q) = 5^ 
and c M (q, t) are the annihilation operators of the particle with momentum hq. By 



substituting Eqs. ( 1122ft and ( 11Q2[) in (1117p . the Hamiltonian of the total system in the 
large-time limit can be recast in the following form 

H = Hq + H int 
Ho = H ete + H F 

4 r 

H ele = J2 d 3 q£q4(q^)c M (q,t), (123) 

11=1 J 
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H F = -.y I ' du I d 3 khojd{(k,uj,t)dx(k,u,t) + b\(k,u,t)b x (k,u,t)} ■■ 

A=l J J 

H int = - J d 3 k(J*(k, t) ■ A(k, t) + h.c.) = -ec J rf 3 x J d 3 k^*a • A(k, t)ip 



- ^ V Y (d 3 k(d\ HdJ g(u) 

txt(q)a-s A (k)v(q-k) 
x — 2 7 x ~ 2 , 2 7 x S(q)cy(q-k)& A (k,a;,0)e -/i.e.). 



(124) 



The unperturbed Hamiltonian if = H eie + ifp has the eigenstate 

| ele + rad) =\ ele) \ rad) 

which are the direct product of the eigenstates of H de and H F . we again apply first order 
perturbation theory to treat the transition probability per unit time for a free Dirac particle 
of momentum hq to emit a photon of momentum ftk and energy hoj thereby changing its 
momentum to %{<\ — k) 

IWk = y |(l k | (q-k| H int | q) | 0) | 2 6(y/h?cW + m 2 c 4 - v 7 ^ 2 |q - k| 2 + m 2 c 4 - hJ) 

(125) 

where the argument of the Dirac 5 function displays the conservation of energy. We may 
proceed to calculate the energy loss per unit time as we did in the previous section but with 
a little modification here. The sum over final states must include a sum over the final spin 
states of the particle with positive energy /i — 1,2, also we average over the initial spin 

states 



dw i :; 2 



dt 2 

A=l fj,,fj,'=l 



CO 



y y d 3 k doohuT^^. (126) 
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In order to calculate the above equation we must evaluate the following sums 

2 2 

S = \j2tl K(q)«-e A (k) V (q-k)| 2 (127) 

A=l n,fj,'=l 



1 2 2 

K(q)c*-s A (k) V (q-k)| 2 . (128) 



A=l n,n'=\ 

For this purpose we introduce the annihilation operators 



and we obtain 



A(q) = ^ TPql (129) 

2 \Eq\ 

A(q-k) = Ca ' (q - k) + /3m f + |£?q - k| , (130) 



5 = ^r[(a-e A (k))A(q-k)(a-e A (k))A(q)] 



1 mV , 2(e A (k)- Vl ) 2 Vl -v 2 

2 1 |£ q | |£ q - k | * ~^~ } 



5' = itr[(a.s A (k))A(q-k)(a.8 A (k))A(q)] 

lr, m 2 c 4 2(s A (k) ■ vQ 2 Vl -v 2 

- 2 jl mm ^ { } 

where we have used v = ftc 2 q/i? q and vi and v 2 are the velocities before an after the 
emission of the photon respectively. The sum over polarizations can be carried out as in Eq. 
(|52|). The result is 

V ?/, 2.x , In 9 / ow, 9~7^7 Vx-V 2 . 



S = S' = -±(1- cos 2 0) + -{1 - J(l - v 2 /c 2 )(l - v 2 /c 2 ) - -L^} (133) 
where again is the angle between q and k given by 



CO . fou ^c 2 , / f 2 . 

We mention here that in the classical theory only the first term occurs on the right-hand side 
of the equation ( 11341) and the second and third terms are a consequence of non relativistic 
and relativistic quantum theory respectively. These terms are small since the wave length 



of the electron is much smaller than the photon's wave length j39|-[40|. Also the second 
term in Eq. (I133P is a small correction to the result we found in the preceding section. We 
neglect this term and the rest of the calculation is similar to what presented in the preceding 
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section. The only difference here is that Eq. f |134j) must be used instead of Eq. (IllOp . The 
result is 



dW e 2 v /' +00 J , ; f + °° , M to 2 M hu ,k 2 c 2 1N I v 2 l2 , 
dkk / dwu(l - — [1 + — - -(— - lWl - -] 2 ) 



(it 27r 2 e Jo ./o t ' 2 ^ 2 2mc 2 to 2 V c 2 " 

x Jm ( 2 / \ i 2 2 7 J • ( 135 ) 

Comparing this expression and Eq. (11 111) with its classical counterpart, Eq. fl52|) . it is seen 
that the only difference is in the argument of the 8 function. In fact, we note that the 
integrations over k and u> in classical equation ( 1521) diverge but its nonrelativistic and rela- 
tivists counterpart that is Eqs. fillip and (11351) have no divergent behavior. To appreciate 
the physical significance of this subject, we consider a transparent and nondispersive mag- 
netodielectric medium. We note that Eqs. (IllOp and (I134p for < v < c and n > 1 provide 
a cutoff in frequency u < w ™ onre/ > re/ which in non- and relativistic regimes are respectively 
given by 

u,r ra = 2mC2{ f, - l) (136) 



and 



, = 2,^(^-1) 

where /3 = v/c. These cutoffs are still bounded above by the electron energy mc 2 or 



mc 2 / yl — /3 2 , respectively. Thus, we can use one of these cutoffs for the range of inte- 
grations to obtain a physically acceptable result for the classical radiation intensity. 
It is easy to show that in a transparent magnetodielectric medium Eq. fl 1 3 5 [) becomes 

(138) 

which tends to the correct relation in the classical and extreme relativistic limits. 



D. Finite temperature Cherenkov radiation in quantum regime 

Our considerations so far have been applied to zero temperature. The inclusion of tern- 

n 

perature may be done in the usual manner [27J . In this case the transition probability (I125P 
for a free Dirac particle of momentum h(\ to emit a photon of momentum Kk and energy hw 
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thereby changing its momentum to h(q — k) is obtained as 

IV q - k = y \H mt \ 2 (iV k +l)(l-n F (q-k))5(^V + m 2 c 4 - V^c^q - k| 2 + m 2 c 4 -^), 

(139) 



where 



^ = ^It _ i ( 14 °) 



and 



n F (q-k) = — 1 . (141) 

eV /fi 2 c 2 |q-k| 2 +m 2 c 4 /fc s T _|_ ^ 

Here the factor iV^ + 1 comes from the phonon creation operator for transition from the 
initial state with 7i k photons to the final state with + 1 photons. In fact, we use the 
thermal average of which is N^. This photon emission process take place as stimulated 
and spontaneous emission. Similarly, the factor 1 — np{(i — k) is the probability that the 
electron state q — k is empty, so that the operator c^(q — k) can create an electron in that 
state. There is also a factor np(q), which is the probability that q is occupied with an 
electron that is unity. But in calculating the total energy radiated as Cherenkov radiation, 
we have to omit some important processes. They arise from other electrons in the system 
with the same spin state. These electrons can not be found in the state q since our electron 
is occupying it already. Thus the other electrons of the system will reduce the transition 
probability (I139p and we must subtract the transition probability due to the presence of 
other electrons from Eq. fl 1 3 9 f) [23]. Therefore by using Eqs. ( I126p - (1135p . the total energy 
radiated in finite temperature is found to be 



dW e 2 v f +co „, f + °° , n u 2 M hu t k 2 c 2 1N / v 2 , 2 , 
dkk / duu(l - -Wl + T^—ji—r - l)\/l - -d 2 ) 



dt 27T 2 e Jq Jq v 2 k 2 2mc 2 cu 2 V c 2 

— u 2 e(u) + k 2 c 2 K>(uj) 



where 



F T (u) = (iV k + l)(l-n F (q-k))-iV k n F (q-k) 

e hui/k B T e \Eq-hui\/k B T _ e -frw/k B T 



(143) 



e huj/k B T _ I !• e \E^-hu\/k B T ' 

The last equation differs from the zero temperature Equation fl 1 3 5 [) only by the multiplicative 
factor F t (uj) which has the following asymptotic behavior at low temperatures k B T « frw 
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and at high temperatures ksT » respectively: 



Ft(co) ~ 1 



F T (u) 



3l 

2hco 



mc 



2hu\ l - (v/cY 



It is easy to show that in a transparent magnetodielectric medium, the finite temperature 
generalization of the equation fl 135[) becomes 

dW e 2 v 



dt 47reoC 2 



r+oo 

„■ JO 



(144) 



which in a nondispersive medium tends to the result of 15]. 



IV. CONCLUSION 

I — I 

In this paper, we have generalized a Lagrangian introduced in [llj to include the external 
charges. In this formalism the medium is modelled with two independent collections of 
vector fields. The classical electrodynamics in the presence of a polarizable and magnetizable 
medium is discussed and the susceptibility functions of the medium are calculated in terms of 
the coupling functions. The energy loss of a point charged particle per unit length emitted 
in the form of radiation, called Cherenkov radiation, is obtained in both zero and finite 
temperature in classical regime. A fully canonical quantization of both electromagnetic field 
and the dynamical variables, modelling the medium, is demonstrated. In Heisenberg picture, 
the constitutive equations of the medium together with the Maxwell equations are obtained 
as the equations of motion of the total system. The wave equation for the vector potential 
is solved. It is shown how vector potential operator in this theory can be expressed in terms 
of the medium operators at an initial time. The consistency of these solutions for the field 
operators are found to depend on the validity of certain velocity sum rules. It is also shown 
how this scheme is related to the damping polarization and phenomenological quantization 
theories. The large-time limit and quantum mechanical perturbation theory is applied to 
treat the finite temperature Cherenkov radiation in the domain of the non-relativistic and 
relativistic quantum regimes. The total energy radiated per unit time is calculated which 
is consistent with both classical and extreme relativistic limits. The approach is based on 
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a Lagrangian formalism and magnetic properties of the medium and the relativistic motion 
of the particle are included. This model can be applied to the case of Cherenkov radiation 
in a nonlinear medium which is under consideration. 



Appendix 

In this appendix we give the analogous expression for the vector potential in the quantum 
domain. Furthermore, we evaluate the time- dependent coefficients for the vector potential 
operator. By using Eqs. ( 155]) and the wave equation ( ETj) . can be written as 

/WL+ \k\ 2 A + fi -^J dt' X e{t-t')A(k,t')-Ho\k\ 2 J dt'x m (t-t')A(k,t') 
dp N± 

= Vo^f- + /^k x + no J x . (A. 1) 

We find the full time dependence of the vector potential by taking the inverse Laplace 
transform. The inverse Laplace transformation as defined before is a contour integration 
over the Bromwich contour. After transforming to frequency variables we obtain 

A(r,t) = J2j rf 3 kW / ^^(r ? (fc,t)a A (k,0)e tk - r + / i .c.)e A (k) 

A— 1 * 



A— 1 V 

+ ^E/ rf3k / duJ^^( 1 (u,k,t)b x (k,0)e^ + h.c^s x (k) 



2 



+ f C(A; ' t)(i " a5(r " ra) + S ( r ~ r ^ ' e ^( k ) e A(k) (A.2) 



2 

a A=l 



where 



- L 12 dae " ^ = - '^-'-"^ 

(A.3) 



2vrz 7„ JOO (s + /u;)(s 2 e(s) + k 2 c 2 k(s)) —uj 2 e(u) + k 2 c 2 n(u) 

(A.4) 
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ICO 



(s + tio)(s 2 e(s) + k 2 c 2 K,(s)) —uj 2 e(uj) + k 2 c 2 n(uj) 

(A.5) 



I St 1 

C(/c, t) = — ds—-—^—-— = —Y j m ( e -*V_X_). (A.6) 

sv ' ; 2m J_ loo s 2 e(s) + k 2 c 2 K(s) kc ^ v c/c(fy) v ; 

The solution of the wave equation is the sum of a transient and a permanent part. The 
latter are expressed solely in terms of the initial medium operators. Long after the initial 
time, vector potential operator will be a function of the medium operators alone since it has 
poles on the imaginary axis in the complex s-plane. With the velocity sum rules discussed 
in Sec. 2, one can see that the coefficient rj(k : t) equals 1 at time t — whereas the other 
coefficients have the initial value 0. 
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